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Abstract 

The semi-inclusive B-meson decay into spin-singlet D-wave 2 h charmonium, B — > w C 2 + X, is 
studied in non-relativistic QCD (NRQCD). Both color-singlet and color-octet contributions are cal- 
culated at next-to- leading order (NLO) in the strong coupling constant a s . The non-perturbative 
long-distance matrix elements are evaluated using operator evolution equations. It is found that the 
^■l color-singlet 1 L>2 contribution is tiny, while the color-octet channels make dominant contributions. 



The estimated branching ratio B(B — > rj C 2 + X) is about 0.41 x 10~ 4 in the Naive Dimensional 
Regularization (NDR) scheme and 1.24 x 10 -4 in the t'Hooft-Veltman (HV) scheme, with renor- 
malization scale [i = mj, = 4.8 GeV. The scheme-sensitivity of these numerical results is due to 

i f8l i [81 

cancelation between Sq and contributions. The /U-dependence curves of NLO branching 

ratios in both schemes are also shown, with // varying from ^ to 2mb and the NRQCD factoriza- 
tion or renormalization scale fi\ taken to be 2m c . Comparison of the estimated branching ratio of 
B —> r} C 2 + X with the observed branching ratio of B — > X(3872) + K may lead to the conclusion 
that X(3872) is unlikely to be the 2 ^ charmonium state r) C 2- 
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I. INTRODUCTION 



One of the missing states in the charmonium family, the f^ 1 -^), is the only missing spin- 
singlet low-lying D-wave charmonium state. Its mass is predicted to be within 3.80 to 3.84 



GeV 



3j, which lies between the DD and the D*D thresholds. The J PC quantum number 
of 7] C 2 is 2 h , thus its decay to DD is forbidden. Therefore, this is a narrow resonance state, 
and its main decay modes are the electromagnetic and hadronic transitions to lower-lying 
S-, P-wave charmonium states and the annihilation decays to light hadrons. Previously, we 
calculated the inclusive light hadronic decay width of the 1 _D 2 state at next-to-leading order 
(NLO) in a s [4j in the framework of non-relativistic QCD (NRQCD). The results show that 
with the total width of r] C 2 estimated to be about 660-810 keV, the branching ratio of the 
electric dipole transition t] C 2 — > ^h c is about (44 — 54)%, which will be useful in searching for 
this missing charmonium state through 7] C 2 — > ^h c followed by h c — > ^ri c and other processes. 

The NRQCD factorization method 5] was adopted in our calculation of r] C 2 light hadronic 
decay. Within this framework, the inclusive decay and production of heavy quarkonium 
can be factorized into two parts, the short-distance coefficients and the long-distance matrix 
elements. A color-octet heavy quark and anti-quark pair annihilated or produced at short 
distances can evolve into a color-singlet heavy quarkonium at long distances via electric or 
magnetic transitions by emitting soft gluons, This color-octet mechanism has been used to 



11 



12 1 charmonium 



remove the infrared (IR) divergences in P-wave [5l-ll0| and D-wave 4j, 
decays. 

Now, we turn to the B-meson non-leptonic decays, which have played an important role 
in discovering new resonances, especially new charmonium and charmonium-like states in 
recent years. The branching fractions of B-meson inclusive decays into S-wave and P-wave 
charmonia, of (9(10 -3 ) to O(10~ 2 )[13], are relatively large. Therefore, we may also expect 
to search for D-wave charmonia in B-meson decays, and in particular to search for the spin- 
singlet D-wave charmonium rj C 2 in B — > r] c2 + X. Like the charmonium light hadronic decay, 
charmonium production in B-meson semi-inclusive decay may also be factorized in NRQCD 



as 



T(B ^H + X) = ^C(6^cc[n]+X)(C H [n]), (1) 

n 

where the sum runs over all contributing Fock states. The short- distance coefficients (7(6—7- 
cc[n] + X) can be perturbatively calculated up to any order in a s ; while the long-distance 



matrix elements (0 H [n]} should be determined non-perturbatively. One may refer to [hJQ] 
for more discussions on the feasibility of Eq. ([T|). 

S-wave and P-wave charmonium production in B-meson semi-inclusive decays have al- 



ready been studied by many authors in the literature |10l . |14H18| . In 1C| |l4( , it was found 
that the experimentally observed branching fractions for J/ip and ip' could be accounted 
for by NLO calculations, while for Xci an d Xc2 the branching ratios were still difficult to 



explain. In |19j, the branching fractions for D-wave charmonium production in B-meson 
semi- inclusive decays were calculated to be of C(10~ 3 ) in NRQCD at leading order (LO), 
where the NRQCD velocity scaling rules were used to estimate the long-distance matrix 
elements. Similar results but somewhat larger branching fractions were also obtained in 
[20! . However, the NLO QCD corrections are found to be very important in many h eavy 



quarkonium production processes, e.g. in e + e~ annihilation[21], hadropro duct ion |22|, |23| . 
and photoproduction|24|]. Moreover, the velocity scaling rules are too rough to give a quan- 
titative estimate for the long-distance matrix elements. Therefore, for D-wave charmonium 



20| . a NLO calculation and a 



production in B-meson semi-inclusive decays, aside from 19 
better estimate for the matrix elements are necessary. 

Another important motivation for carrying out this study concerns the long-standing 
puzzle of the nature of X(3872). Previous studies assumed that the quantum numbers 
of the X(3872) were J PC = 1 ++ , and this was supported by a number of measurements. 
However, the new BABAR measurement of X(3872) — > J/^7r + 7r^7r 25( favors the negative- 



parity assignment 2 + . Nevertheless, people stil 



X(3872) strongly disfavor the 2 + assignment 26 



argue that the observed properties of 
Recently, [30[ proposed that the 



-2i 



angular distributions of decay products could be used to distinguish between the 1 ++ and 
2 h assignments of X(3872). In this paper, we will further clarify this problem by calculating 
the l D 2 charmonium production rate in B-meson semi-inclusive decay. We will compare the 
calculated branching ratio B — > i] C 2 + X, with the experimental measurement of Br(B — > 
X(3872)K), and then discuss if X(3872) can be the 2 _+ charmonium rj C 2. 

The paper is organized as follows. In Sec. Hi] and 1111} decay widths of four contributing 
Fock states at tree and one-loop levels are calculated both in QCD and NRQCD, and finite 
short-distance coefficients C(b — >■ cc[n] + X) for different components cc[n] are obtained 
respectively after matching between QCD and NRQCD. Computation methods adopted 
in real and virtual corrections are discussed too. The long-distance matrix elements are 
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estimated using operator evolution equations. In Sec.lIVt numerical results are given and 
analyzed. And finally the possibility of assigning the r] c2 as X(3872) is discussed. 



II. LEADING-ORDER (LO) CONTRIBUTION 



We use the same description as in [lo[ Q. The weak decay b — )■ cc + s/d occurs at energy 
scales much lower than the W boson mass mw Integrating out the hard scale and making 
Fierz transformation, we finally arrive at the effective Hamiltonian 

6 



y/2 



cq 



q=s,d 



r q 1] ( M )e» 1 ( M ) + q 8] ( M )c» 8 ( M ) 



i=3 



where the cc pair is either in a color singlet or a color octet configuration, denoted by 0\ 
and Og respectively, 



1 = [c Tm (1 - 75 )c][&7 M (l - 75)9], 

8 = [cT a ^(l- l5 )c][bT a Y(l-l 5 )q]- 



(3) 



(9 3 _ 6 are the QCD penguin operators [3l|. C[i](/x) and C[ 8 ](/i) are the Wilson coefficients of 
0\ and 08, and related to another group of coefficients C + (/i) and C_(/x) through 



= 2C + (/i)-C_(/i), 



(4) 



At LO, expressions for C±(fi) are 

with the one-loop anomalous dimension 

,(0) 



(raw) 



(5) 



7i u; = ±2(3=Fl), 



and a. s 



47T 



/3 ln[/iV(A^)2 



(6) 



(7) 



where f3 = 11 - §JV). We choose m w = 80.399 GeV[13|, m z = 91.1876 GeV, m b = 4.8 GeV, 



Nf = 4, and Aq^. d = 128 MeV for four flavors to adjust a s (rriz) to be 0.119 for five flavors. 
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Only four configurations contribute to r] c2 production at LO in v, the velocity of heavy 
quark or anti-quark in charmonium rest frame: 

\ Vc2 ) = 0{l)\ l D [ 2 1] ) + Oiv^P^g) + O^S^gg) + ■■■. (8) 

With the Fock state expansion Eq. (jSJ), we have 

r(6 -> r] c2 x) = r(6 + r(6 ->• l sfx) + r(6 -> ^x) + r(& -> ^'x) 



-5 r^v ^1 I J" 



(Oi^Sq)), (Os^Sq)), (Os^Pt)) and {Oi( 1 D 2 )) are the production matrix elements of four- 



fermion operators defined in 



32|: 



O 1 ( 1 S ) = xH{^ H a H )^X, 
OsCSo) = ^(oW)^; 

OiCiPa) = (fl^an) V^'x, (10) 

where = B - S and 5« = (-i) 2 (^*^i - §S 2 ^). 

We use Wolfram Mathematica 7.0.1.0, feynarts-3.4, and FeynCalc 6.0. At tree-level, 
the coupling vertex structure C7 M (1 — 75)0 restricts possible J PC numbers of charmonium 
states. Matching amplitudes in both QCD and NRQCD at LO leads to finite short-distance 
coefficients 

CfS®) = r Cft3(l - rjf, 



ces®) = r c^-(i- v T, 



,« 9 

1 - 1 0<^'[8 

cfpf 1 ) = 0, 

CCD?) = 0, (11) 



where 



r ° " 2167r(2m c ) ' 17 (12) 
and |K S | 2 + \ V c d\ 2 ~ 1 has been used. For the LO Feynman diagram, see Fig. [JJ. The strong 
dependence on renormalization scale \i of 8 , (//) at LO causes the results in Eq. (ITT]) 



(9) 
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(tree) 



s/d 



FIG. 1: LO Feynman diagram of b —> cc + X. 
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FIG. 2: LO //-dependence curves of Chg\([jt,). The solid line denotes Cm(/z), and the dashed line 
C[ 8 ](/i). Ratios of g j(/i)/C^ 8 ](m&) as functions of \i are also shown. 



unreliable (see Fig. [2]) and calls for higher order corrections. The QCD Penguin operators 
in Eq. ([2]) also contribute to non-zero tree-level decay width, although their contribution is 
tiny due to the smallness of C$-§{pi). We will neglect their /^-dependence and adopt those 
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14] . for they chose the same values for m&, ^qcd as ours. C 3 (m b ) = 0.010, 
0.007 and C 6 (m b ) = -0.028. Together with Ch?(m b ) = 0.42 



values given in 
C A (m b ) = -0.024, C 5 (m b 
and C^°(m 6 ) = 2.19, the Pen euin contribution is 

3(C 3 -C 5 ) + C 4 



5 P [ l S [ ' ] ] 



°[i] 



0.06, 

c 4 — c 6 

0.007, 



(13) 



6 



which add corrections to tree-level short-distance coefficients in Eq. {TT 

<2 
'[1] 



cfpf 1 ) = o, 

CCD?) = 0. (14) 

III. NLO CALCULATION AND DIVERGENCE CANCELLATION 
A. Real Corrections 

Gluon mass regularization is adopted in our calculation, therefore 75 matrix can be treated 
in 4-dimension. Real correction figures are in Fig. [3]. Divergences are separated from the 




V 'V 
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(n) (r 2 ) 
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b , . s/d b s/d 

(r 3 ) (n) 



FIG. 3: Real correction Feynman diagrams of b — > cc + X. 



finite parts in the amplitude squared. Two kinds of divergences appear: the soft and the 
collinear. Three divergent regions exist: soft, soft-collinear and hard-collinear. Take 1 5 , q 1 ' for 
example. In the soft region, the gluon connected to the incoming bottom quark turns soft, 
i.e., its momentum goes to zero {(r\) of Fig. [3]); in the soft-collinear region, b-quark gluon 
turns soft and at the same time s/d-quark gluon is collinear with the outgoing s/d quark, or 
their momenta are parallel to each other ((ri) and (r 2 ) of Fig. [3]); and in the hard-collinear 

7 



region, s/d-quark gluon runs parallel to the s/d-quark ((r 2 ) of Fig. [3]). IR divergences in 
(7-3) and (r^) of Fig. [3] cancel each other. We take the following parametrization 

6(pi) ->■ c(p 4 ) + c(p 3 ) + s/d(p 5 ) + </(pe), (15) 

and the quark propagators in four quark lines have denominators 

Nt = -2 Pl -p 6 +pl, 
N 4 = 2pi -pe+pl, 
N 3 = 2p 3 -pe+pl, 

iV 5 = 2p 5 -p 6 +pl J (16) 
respectively. For 1 S'q 1 ', ^3 = p^ and N 3 = N 4 . Divergent terms are extracted before doing 



phase space integration: 

soft terms 
soft-collinear terms 



1 

ivf' 
1 



NxN< 



5 

1 1 



hard-collinear terms : ~ — , ~ j-^. (17) 

Some of the hard-collinear terms are seemingly divergent but finally contribute to the finite 
parts. The Mandelstam variables are 

s = ( Pl -p 6 ) 2 , 
t = (Ps + Pe) 2 , 

u = (Pi-p 5 ) 2 , (18) 

and 

u = Am 2 + ml + X 2 — s — t, (19) 
with A the non-zero gluon mass. Rescaling all the dimensional variables with respect to mj 

m c = -yV^ 
A = m b ^, 

(20) 



S 



and 



s = ml(l-y + £), 

t = ml(l-x + T]), (21) 

we finally arrive at the amplitude squared expressed using dimensionless variables x, y 
instead of s and t. Upper and lower limits of x and y are derived from those of s and t via 
Eq. (J2I]) 

y max = 1 + C - 777-^ r(2r? -x+ v^ 2 - 4r/)(-2 + 2£ + x - v^ 2 - H , 

4(1 + rj — x) 

y mi n = 1 + £ - -uTV 1 V ( 2r ? ~ x ~ V^ 2 -4r/)(-2 + 2£ + x + ^x 2 - Ar]) , 

4(1 + T] — X) 

•Ernax 1 ^ ~\~ V i 

Phase space integration over x is a little bit complicated, and the Euler transformation is 
needed by introducing a new integration variable 



tt = J X -^4 (23) 



to replace x and its integration limits 



Umax 



'77-2^7-e + l 



r/ + 2^-e+l 

ttmin = 0. (24) 

Divergences in (r 3 ) and (r 4 ) of Fig. [3] can not cancel each other for ^q 8 ', which makes diver- 
gent terms more complicated. They also produce the only IR pole, the residual divergence in 

lp[8] 

, which can be cancelled by absorption into the redefinitions of non-perturbative matrix 

elements of and 1 S'q 8 ' states. Furthermore, there is no divergence in real correction of 
i D [i] 



B. Virtual Corrections 

In virtual corrections, IR divergences, soft and collinear, are regulated with non-zero gluon 
mass like in real corrections. Ultraviolet (UV) divergences are dimensionally regulated at 
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the amplitude level before projecting the free charm quark pair onto certain charmonium 
bound state of particular angular momentum and color. Virtual correction figures are in 
Fig. [3J. Each diagram in Fig. [1] has an loop integration over gluon momentum q. For 





to.) 



s/d 




s/d 




(«s) 



s/d 



7 



FIG. 4: Virtual correction Feynman diagrams of 6 —)• cc + V. 



example, in (t>i) the UV divergent loop integration has the form 



d D q 



q PqP 



(25) 



(2n)D (q 2 - A^)(( Pl - g)2 - mg)((p 4 - g)* - m 2 ) ' 
and the UV divergent term comes only from the region when q — > oo 

f d D q qPqP' 

J (2vr) D q 2 ■ q 2 • g 2 ' 1 j 

which is proportional to the D-dimensional metric tensor g pp ' . Thus corresponding fermion 
chain in (vi) reduces into 

r M 7p7 Q ®TVr M . (27) 

is the short form for electro- weak vertex 7 M (1 — 75). UV divergent term extractions from 
structures like above are carried out upon using the Fierz transformations 

7p7a^®7 P 7 Q r M = (W + AXeuv^^T^ + Ex, 
r M 7p7 a ®7°7 P r At = (4 + 4Ye uv )r,®r p + E Y , 

® 7p7 Q r^ 7 °7 P = (4 + 4Ze uv )T, ® + E z , (28) 
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where the scheme dependence of 75 is fully extracted and contained in scheme- dependent 
variables X, Y and Z, 



NDR scheme : X = -1, Y = Z = —2; 
HV scheme : X = -1, Y = Z = . (29) 

Hence, the 75 matrix in can still be kept in 4-dimension when evaluating the trace 
formalism. Evanescent operators Ex, Ey and Ez exist only in D 7^ 4 dimensions but 



vanish in D = 4 31] . Therefore they make no contribution to the decay widths, and can 
be discarded throughout the calculations. Again for the 1 5'g 1 ^, self-energy diagrams of (v 3 ) 
and (vq) can only exist for color-singlet electro-weak vertex, i.e., only CW/i) appears. On 
the contrary, the other four diagrams (^1,2) and (^4.5) can only have CW/i) electro-weak 
vertex. Those six diagrams only couple to the tree diagram with C[i] (/z) vertex, contributing 
to C|jj(/i) and C[i](/x)C[ 8 ](/x) terms, respectively. IR divergence of (%) cancels that of (v 2 ), 
and (1*4) cancels (t> 5 ). 

Adding self-energy diagrams in Fig. [5], one can remove UV divergences in (t^) and (vq). 
Explicitly, 



s/d fa s/d 

(-0 M 



V 




s/d 6 s/d 

(*0 W 



FIG. 5: Self-energy correction Feynman diagrams of b —> cc + X. 



(f 3 ) + (si) + (s 2 ) = UV finite, 
(ve) + (S3) + (s 4 ) = UV finite, 



(30) 
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where 



(si) = --i (4na s )N t (m b ) 
4 

(s2) = --i{^a & )N t (m b ) 
4 

( s 3) = --^ (4vra s )A^ e (m 6 ) 



(s 4 ) = --i(4ira s )N e (m b ) 



2e uv 

^ + 3 log (|)_ logK) _ 2 



1 
1 



(tree), 
(tree), 



- log(0 - 2 



, MO 1 

26^ 2 4 



(tree), 
(tree), 



(31) 



with iV e (m 6 ) = i (47r) e ^ 2 r(er/y + 1) ^ No virtual corrections to ^f 1 and 1 D [ 2 1] 

exist accurate to NLO in a S) because of their vanishing tree-level amplitudes. This leads to 
a convenience directly that computation is reduced significantly, (vi) + (r> 2 ) + (1*4) + (1*5) is still 
UV divergent, which needs operator renormalization, i.e., to subtract the term proportional 



to — 7b + ln(47r) or equivalently make the replacement 



1 



— > — ln(47r). 



(32) 



7b is the Euler constant. To summarize our renormalization procedures. First, make mass 
renormalization for charm, anti-charm and bottom quarks — > m = + m ct (No such 
operation is needed for strange or down quarks which are taken to be massless in this paper.), 



m ct = -i{4ira s )N e (m b ) 



+ 4 



(33) 



second, add the self-energy diagrams of external quark lines; finally, do operator renormal- 
ization explained above. 



C. Residual Divergence Cancellation 

We then demonstrate how the residual IR divergence is cancelled. At NLO in a s , on the 
QCD side, 

F(b ->■ 1] c2 X) = C'( 1 5'^ 1 )^^ e+c , ouZom6 (Ci( 1 5'o))B or -„ + SQ ] )®^ e+Coulomb (0 8 C S )) Born 

, ^yl p[8h QCD (^8 Q-Pl)) Born , n( \ n [lKQCD Born /nA\ 

+ )soft — 2 + U l U 2 ) finite ~l > l d4 J 



while on the NRQCD side, 



m" c 1 ml 



T(b^ Vc2 X) = CCS [ 1] ) NR (OiCSo)) NR + CCS [ 8] ) NR (OsCSo)) NR 

+ c( 1 p1 8] ) nr ^^^^- + c^ptj ] ) NR ^ 1 ( 1 -° 2 )) JVR 



mi 
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The subscript Coulomb or soft means having Coulomb or soft pole. NRQCD operator 
mixing of 1 5'q 1 ' 8 ' and 1 p] 8 ' is shown in Fig. [6]. Similar for 1 P 1 ' 8 ^ mixing with l D l 2 ] . And the 



loH.BJ 




lc [l,8] ioP-,8] 



(mi) 




lc[1.8] 









^1 


Annnnnnnnnnra 






o^jooooooooooX 


{m 3 ) 


, ls [i. 8 ] 


ls [i. 8 ] , 


(m 4 ) 



lc[l,8 



lt?[l,8] 1d|1,8] 
°0 °0 



1<j[1.8] 




(m 5 ) 



lt?[1.8] lcj[1.8] 
"0 "0 



1q[1,8] lc[l,8] 



(me) 



FIG. 6: NRQCD operator mixing of 1 S [ q' S] and 1 p[ 8] . 



non-perturbative matrix elements up to NLO in a s are 

(OiCs )) NR 

(O 8 CS )) NR 



(Oi^So^Born + (0\ C Sq) ) Coulomb ~ T^O 11 ~ + o)(~o~ 



^in A 2 l v 16 (0 8 ( 1 P 1 )) flo rn 

4tt I %i + 3 n 3 J m 2 



Born Coulomb 
a A 2 1 16 / <Ol( 1 Pl)>Bom , p (Q 8 ( 1 i 3 l))Bom 

4^ (ln ^ + 3 )( y } ^ 2^ + Bf ml 



(0 8 CPl)) NR = (0 8 CPl))Born + (^(^l)) 



Coulomb 



a 



~ 7z(ln-T+o)( : 



A 2 , l w 16. / (0 1 ( 1 D 2 )) Born , (0 8 ( 1 J D 2 )) Born 



4vr v >| ' 3 yv 3 ' 2N c m 2 c 



C 



mi 



(35) 



Bp = j^. The Coulomb singularity in (m.5) and (wig) of Fig. [6] is extracted and related to 
the tree-level matrix element in the following way 



(0[ n ](cc)) Coulomb = C[ n ]^- (0[ n ](cc)} B orn, 



(36) 
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with the color factor 



Cp — 3, n=l color-singlet cc 



C W = A (37) 



•^r = — ^ n =8 color-octet cc. 



leading to 



T{b^ Vc2 X) = CCS^iOlCSo^Born + CCS^Bor^CSo)) 
- C{ b J Bom T- (m -y + -)( — 



Coulomb 



An y ~" /i\ 3 3 m 2 

Born Born Coulomb 

r( lq[8U a sn A 2 1^16 (0 8 ( 1 Pi))Bom 



'4ji 3 3 m 2 . 

I ^^lplS^NR i^S^Pl^Born | ^ljjM jiVfi (^1 )fiorn 



Matching Eq. ( )34|) and Eq. (1381) . one can get the finite short-distance coefficients accurate 
to one-loop level 

r^nq^hNR _ r (iq[t]\QCD 

U{ d Q ) — U{ D Q ) finite , 

rvi c^hW-R _ rvi q[8KQCD 

l D J — °( D ) finite ' 

+ C( 1 4 1] W n ^(ln4 + ^)(T) 



^Wn^Cln^ + ^y)^, 



C(^W)^ = CCD?)fJ e . (39) 

Coulomb singularities in C{ l Sl ] ) QCD and C^ 81 )^ and soft divergence in C( l Pf ] ) QCD 
are absorbed into the long-distance matrix elements (Oi( 1 Sq)) nr and ((9 8 ( 1 5'o)) Afi? . There 
is no residual soft divergence in real correction to because of the absence of tree-level 
amplitude of 1 P 1 . Considering its vanishing virtual correction, the NLO correction to 1 D^ i 
is finite. One-loop level short-distance coefficient can be expressed in the common form 



C(b -> cc[n] +x) = T ^-(C( 1]gi [n] + 2C [1] C [8] g 2 [n] + C( 8] g 3 [n\), (40) 
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and gi[n], and g 3 [n] of^iVsi 81 and ^Pj 81 were calculated in |l(J, ll4J ■ We list them in 
Appendix. B. For 1 D 2 1 \ our results are new: 

9i[ l Df] = 0, 
92[ l Df\ = 0, 

9z[ l Df\ = ^ {2 V 3 - 9v 2 + 18?7 - 61ogfa) - ll) . (41) 
D. Evaluation of long-distance matrix elements 

Due to lack of experimental information on the matrix elements of D-wave operators, 
we can not extract them from experiments and have to invoke some theoretical estimates. 
The color-singlet matrix element (Oi( 1 D 2 )) may be determined by potential models with 
input parameters, while the color-octet matrix elements may be estimated using the oper- 
ator evolution equations. Matrix elements (0 8 ( 1 Pi)) NR , (Oi( 1 S )) NR and (O s ( 1 S )) NR are 
renormalized in NRQCD, and thus have /i a- dependence, and this can be explicitly shown 
by deriving the quantities on both sides of Eq. fl35l) with respect to fi^: 

diO^So))™ a s 32(O s ( 1 P 1 )) Born 



2 



<iln/iA 4-7T 3 m; 

d((9 8 (%)>^ a s 32 R (0 8 ( 1 P 1 ))b 

-±s F - 



2 



ciln/iA 47r 3 m 

d(O s ( 1 P 1 )) NR 0*32 (Q 1 ( 1 D 2 )) Born 

rfln/^A 4tt 3 F 2N c m 2 c ' 1 ' 



Eq. fj4*2|) has the same form as Eq. (45) in (4], where the IR divergence is regularized in 
dimensional regularization scheme. This is because the operator evolution equations have 
nothing to do with the IR divergent parts. The solutions are 

where we take m c = 1.5 GeV, b = ^ - ^, C A = 3, iV>=3, A^g D = 153 MeV for LO, and 
Aqcd = 399 MeV for NLO. 

The initial matrix elements like (Os{ l Pi){fJ>A ))a.t starting scale fi\ = m c v, where 
v 2 = 0.25, are eliminated. One could refer to |4j for reasonability of doing so. The 
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evolution equation method for determining the long-distance matrix elements has been 



used in estimating the D-wave charmonium state light hadronic decay width and h c decay 
width . |33|. For h c , the evolution equation could give a prediction for light hadronic 



decay width within about 30% error when compared to experimental extraction [33J. That 
means the operator evolution equation is a good method to evaluate the P-wave long-distance 
matrix element, and can be extended to D-wave case, which is lack of experimental data. 



IV. RESULTS AND DISCUSSIONS 

The long-distance CS D-wave matrix element is related to the second derivative of the 
radial wave function at the origin 

151/?" (Q)\ 2 

(O l {n 1 D2)) = (2J+l)(n 1 D 2 \0 1 (n 1 D 2 )\n 1 D 2 ) = 5 (2N C ) 1 f { )l , (44) 



where N c = 3 and B-T potential model input parameter \R'{ D (0)\ 2 = 0.015 GeV 7 [34( for 
charmonium. Before giving the final results, we have to first deal with the NLO Wilson 



coefficients C[i](/x) and C[g\(n). The expressions for C±(/jl) up to NLO in a s are given in [35] 



c±(aO 

with 



, „c°> 

a s (M w 



a s (n) \ Air J \ Att 



(o) a (i) 
7± Pi 7± 



2/3 2 2/3 ' 

B ± = ^(±11 + k ± ), (46) 



and the one-loop and two-loop anomalous dimensions 



7 i 0) = ±2(3=Fl), 

7 « = ^(-21 ±\N } -2^ ± ). (47) 



The scheme-dependent k± are 



0, NDR scheme, 
«± = <( (48) 
=F4, HV scheme. 
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TABLE I: LO (Eq. CE|)) and NLO (both Eq. ([11} and Eq. flU])) short-distance coeffecients of four 
subprocesses, with Tq removed. Results for both NDR and HV schemes are listed. The QCD 
renormalization scale [i takes values from ^ to 2m;,, where m b = 4.8 GeV, m c = 1.5 GeV. 



Fock state 


LO 


NLO NDR scheme 


NLO HV scheme 


A* 


m b /2 


m b 


2m b 


m b /2 


m b 


2m b 


m b /2 


m b 


2m b 


i D m 











0.0028 


0.0020 


0.0015 


0.0026 


0.0018 


0.0014 


lp[8] 











-2.058 


-1.545 


-1.289 


-1.880 


-1.390 


-1.150 




0.0458 


0.2130 


0.4330 


-0.2102 


-0.3978 


-0.4892 


-0.0633 


-0.1950 


-0.2629 


lap] 
°0 


8.803 


8.065 


7.566 


12.856 


11.217 


10.169 


13.490 


11.529 


10.287 



Note here an additional factor — y should be included in B± in the HV scheme. ft and ft 
are in the NLO expression for a s 



4tt 



ft ln[m(/, 2 /A^)] 
/3 2 ln(^/A^ D ) 



(49) 



with Aqcd = 345 MeV, ft = 11 - f JV>, and ft = 102 - f N f . 

LO and NLO short-distance contributions are given in Table [H It is easy to see that at 
renormalization scale fi = m b , the short-distance coefficients in NDR and HV schemes differ 
slightly for the dominant components 1 P 1 ' 8 ^ and 1 Sj^. The long-distance matrix elements 
take the following values 



WD*)) 



rrr. 



mi 



(OiCSo)) 
(OsCSo)) 



0.053 GeV 3 , 

0.0092 GeV 3 , 
0.0036 GeV 3 , 
0.0015 GeV 3 , 



(50) 



(Oaf 1 Pi)) 



where m c = 1.5 GeV and [i\ = m c v = 750 MeV. The long-distance matrix elements 
(Oi^So)) and (Os^Sq)) are sensitive to charm quark mass m c and initial scale /xa - Multi- 
plying the short- distance coefficients shown in Table |T] by the matrix elements in Eq. fl50|) , 
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we get the B- meson semi- inclusive decay width into 77^. Then we can estimate its branching 
ratio using B-meson inclusive semi-leptonic decay rate. That has the benefit of eliminatin. 
the Vfc dependence and reducing the dependence, as was performed in 0, Q, Q 
The theoretical prediction for the inclusive semi-leptonic decay width can be expressed as 



SL 



G 2 F \V hc \ 2 ml 

1927T 3 



[1 - 8z 2 + 8z 6 



24^ 4 log(^)]r ?1 (z) 



(51) 



where z 



37] 



The factor r/i(z), including NLO QCD correction, has the approximate form 



?7i0) 



3tt 



^ + (-y + - 2 )(l--) 2 



(52) 



Using the calculated B-meson semi-inclusive decay width given in Eq. (I5TI) . and the experi- 
mental semi-leptonic branching ratio Btsl — 10.74% [13J, and taking m c and fi\ in regions 
(1.4, 1.6)GeV and (700, 800)MeV, respectively, we finally arrive at the QCD renormalization 
scale /i-dependence curves in Fig. [7] for the branching ratio Br[B — > ri&X] of B-meson 
semi- inclusive decay into 77,^. Note that varying /ia only changes the relative ratios among 
long-distance matrix elements, while varying m c affects not only the long-distance matrix 
elements but also the short-distance coefficients. 




x 

(S 

t 
CO 



-I 



I'---- 



5 6 7 
H (GeV) 



FIG. 7: QCD renormalization scale /i-dependence of Br[B -4 ri&X] in NDR scheme (left) and 
HV scheme (right). The long-distance matrix elements are estimated using operator evolution 
equations, [i ranges from ^ to 2rrib. The shaded zone is for the values of Br[B — > rj&X]. Upper 
bound for solid curves correspond to m c = 1.4 GeV and /i\ = 700 MeV, dashed lines to m c = 
1.5 GeV and /iA = 750 MeV, and lower bound solid curves to m c = 1.6 GeV and /_iA = 800 MeV, 
respectively. 
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When \i is taken to be m& = 4.8 GeV, 

Br(B ^ Vc2 X) NDR = (0.41+^) x 10~ 4 , 
Br(B ^ Vc2 X) HV = (1.24^-23) x i " 4 , (53) 

where the central values correspond to m c = 1.5 GeV and n\ = 750 Me V, upper bounds to 
m c = 1.4 GeV and /xa = 700 MeV, and lower bounds to m c = 1.6 GeV and /xa — 800 MeV, 
respectively. Since the color-octet Wilson coefficient C[s\(n) is much larger than the color- 
singlet one C[i](/i) 

~ 15 , (54) 



the LO decay width is dominated by that of x Sq \ which is proportional to CL{^i). For NLO, 
decay widths of and l D 2 are negligible, and those of 1 P| 8 ' and 1 5 , g 8 ^ are of the same 
order and make most contribution to the branching ratio in Eq. (153]) . but unluckily they 
largely cancel each other. This cancellation is related to our estimates for the long-distance 
matrix elements in Eq. ( 1501) . If without this cancellation, the l sf^ Fock state could give the 
following central values 

Br(B -> 1 S [ Q ] X) NDR = 5.30 x 10" 4 , 
Br(B -> 1 Sf ] X) HV = 5.45 x 10~ 4 , (55) 

which might be regarded as the upper bound of the branching ratio for this process. Fur- 
thermore, we may consider the following uncertainty in the predictions of the branching 
ratio. Since 

C[1] '^ s , (56) 



°[8] 

we might carry out a double expansion in both a s and C[i]/C[ 8 ] simultaneously [15] . In this 
new expansion, terms of different orders scale as follows: 

LO: C[|] ; 
NLO: a s C 2 8] , C[i]C[ 8 ] ; 
N 2 LO: a 2 C[ 2 8 ], a s C[i]C[ 8 ], Cj 2 ] ; 

N 3 LO: ttsC[ 2 8], a 2 C[i]C[ 8 ], atsC^, • • • . (57) 
C 2 fil scales as LO, and c^C 2 , as NLO. a?C 2 i scales the same order as a s C^C[g\ and C 2 j, thus 



should also be considered. Authors of Il_5| did not calculate all a 2 C|] terms, but estimated 
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their contribution by adding a correction term of the same order. The same method with a 



minor modification was adopted in 



10 



, Jjj]. Unluckily, their method can only be applied to 
the color-singlet channels that have non-vanishing LO decay widths, and fails in our case. 
In [l8] the a^C^] virtual contribution from squared one- loop amplitudes was calculated, 
but the real correction was neglected by arguing that the real contribution was phase-space 
suppressed. However, the IR divergent real corrections can not be omitted, as pointed out 
in floL fl^j . Hence, a complete calculation at NNLO in a s might be needed to obtain the 



afCjlj contribution, but this is already beyond the scope of our calculation in this paper. It 
will be interesting to see if the large cancellation of 1 P 1 ' 8 ^ and ^q 8 ' could be weakened after 
including the o^CL contribution. 

We now discuss the possible relation between the semi-inclusive decay branching ratio 
B — > rjtfX and the exclusive decay branching ratio B — > rj^K. Obviously, the latter must 
be much smaller than the former, since the X includes many hadronic states other than the 
kaon. In particular, in the case of B — > rj&X, the dominant contribution comes from the 
color-octet cc channels, which subsequently evolve into r] C 2 by emitting soft gluons which 
then turn into light hadrons such as pions. Whereas the exclusive process B — > rj^K 
requires the soft gluons be reabsorbed by the strange quark in b — > cc + s. This probability 
is apparently very small. As a conservative estimate, we believe the branching ratio of 
B — > i] c2 K should be smaller than that of B — > rj^X by at least an order of magnitude. 
The suppression of exclusive decay relative to inclusive decay is supported by many other 
charmonium states. E.g., the branching ratio of B — > J/ipX is (7.8 ± 0.4) x 10 -3 13|], while 
Br(B+ J/ijK + ) = (1.014 ±0.034) x 10~ 3 and Br(B° J/ipK°) = (8.71 ±0.32) x 10" 4 . 
For Xd, Br{B X dX) = (3.22 ± 0.25) x 10" 3 , Br{B+ Xc\K + ) = (4.6 ± 0.4) x 10~ 4 
and Br(B° — > xaK ) = (3.90 ± 0.33) x 10~ 4 . Evidently, the observed inclusive branching 
ratios are about 10 times larger than the corresponding exclusive one. For Xc2, which 
is similar to r] C 2 because in both cases at LO the color-singlet cc Fock states make no 
contributions, Br(B X c2X) = (1.65 ± 0.31) x 10~ 3 , Br{B+ -»• Xc2^ + ) < 1-8 x 10" 5 and 
Br(B° — > Xc2K°) < 2.6 x 10~ 5 , the suppression of exclusive decay is almost by two-order 
of magnitude. Therefore, we may have a general observation that for a charmonium state 
produced in B-meson decays, the suppression factor of exclusive production branching ratio 
relative to inclusive one should not be larger than 1/10 (including the factorizable and non- 
factorizable exclusive processes). This means BriB — > rj&K) should be at most (9(10~ 5 ), 
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based on our calculation. 

In contrast, for X(3872) the observed branching ratio Br(B — > X(3872)K) x 
5r(X(3872) -> D D n ) = (1.2 ± 0.4) x 10" 4 [13]. Considering that there exist many 
decay modes of X(3872) other than X(3872) -»■ D°D 7i°, we may conclude that Br(B -> 
X(3872)K) is at least 10 times larger than Br(B — > rj^K). Therefore, X(3872) is unlikely 
to be the J PC = 2~ + charmonium state 7] C 2- In fact, for X(3872) the J PC = 1 ++ assign- 



ments of the D°D*° molecule 
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40] are preferred 



or a charmonium- .D /}* mixed state 
jy^many authors, instead of a J PC = 2~ + state (for more discussions see a recent review 



ill)- 



V. CONCLUSIONS 



In this paper, we calculate the semi-inclusive decay width and branching ratio of B — >■ 
T] C 2X at NLO in a s in NRQCD factorization framework. The finite short- distance coefficients 
are obtained by matching QCD and NRQCD, and the non-perturbative long-distance matrix 
elements are evaluated by using the operator evolution equations. We find that at tree-level, 
only the S-wave Fock states 1 5'q 1 ' 8 ' contribute, and the LO decay width is dominated by that 
of 1 S , q 8 ', because of the largeness of the color-octet Wilson coefficient squared CLi^i) over the 
color-singlet one C^(fi). Unlike rj c2 light hadronic decay, in this process, there is no residual 
divergence at NLO of the l D 2 1 ^ Fock state, due to the vanishing tree-level contribution of 
1 p| 8 '. At NLO in a s , 1 p| 8 ' and dominate. Unfortunately, they largely cancel each 
other. This cancellation depends on our method for estimating the long-distance matrix 
elements. As a result, we obtain the branching ratio Br(B — > ri^X) = ( 0.41^Q5g ) x 10~ 4 
in the NDR scheme and ( 1.24^QgQ ) x 10~ 4 in the HV scheme, at fj, — mj. The central 
values correspond to m c = 1.5 GeV and fi\ = 750 MeV, upper bounds to m c = 1.4 GeV and 
/Mo = 700 MeV, and lower bounds to m c = 1.6 GeV and /iA = 800 MeV, respectively. If the 
large cancellation does not exist, the 1 S , q 8 ' could give Br(B — > 1 Sq^X)ndr = 5.30 x 10~ 4 
and Br(B ^^X) 

hv = 5.45 x 10 4 , which could be regarded as the upper bound of the 
branching ratio of this process. The /^-dependence curves of NLO branching ratios in the 
two schemes are also shown, where /z varies from ^ to 2m;, and fi\ = 2m c . Furthermore, we 
estimate the exclusive decay branching ratio of B — > ri&K by considering the suppression 
ratios of exclusive decays relative to inclusive ones for other factorizable and non-factorizable 



21 



exclusive charmonium production processes, and conclude that X(3872) is unlikely to be a 
2~ + charmonium state. We hope that our results will be useful in finding the missing 
charmonium state r] C 2 in experiments, and in further studying 7/ c2 production in B-meson 
exclusive decays. 
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VII. APPENDIX 



A. Covariant projector method. 

In our calculation of short- distance coefficients, the covariant projector method is 



adopted 



421 ] . For any spin-singlet charmonium production in 4-dimension, the covariant 



projector is 

Po,o[P, k ) - ~ 7= / ^=7 — 77 1 > looj 

2V2 A m 7 ~ M Im 



2 i m c y y + m c 

where momentum of charmonium bound state P = p<± + P3- Relative momentum between 
charm quark and anti-charm quark satisfies 

P , 

Pi = — + k, 

Ps = ~ - k. (59) 

Bound state mass M = 2m c , which holds in QCD radiative correction calculations, for the 
relativistic effects are neglected. For more details, one could refer to related contents in (J). 
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B. One-loop level short-distance coefficients of 1 S'q 1 ', 1 S'q 8 ^ and 1 pJ 8 ' Fock states. 

For l Si}\ 

g^Sfp] = -4(1 -7 7 ) 2 ^8Li 2 ( ?7 )-4Z + 41og(l - 77) log^) + ^) 



8(2 5t7)(1 77) 2 log(l 7]) 



+ 20(1 - J]) 2 + ^ ^ ^ ^ - 1677(1 - 77) logfa), 

92 [^] = 4(1 - „)> (Slog (4) " X + Y)- 2dV- 53, + 34)(l-,) 



+ 4t/ 2 log(7/) + 



/i 2 y y 2-77 

8(3-t/)(1 -7/) 3 log(l -7/) 



= 1 (-(1 - 77) (2t7 2 - 77/ + 11) - 6108(7^)) ; (60) 

for 

^il^i 81 ] = -^(l-?/)(277 2 -777 + ll)-81og(77), 
92[S ] = 3(1-77) (31og^-j-X + Yj — 

+ 377 2 log(77) + 6 ( 3 -^(;^ 3 ) ^ 1 -^ , 
^ ?(l-„)-(-41og@ )+f + ^-f -31og 2 (2-77) 

1 — 77 



34) (I-77) 
77) 



+61og(l-77)log(2-77)-61og(2)J +3(1- //) ( 9(// l)Li 



2-77 



18Li 2 ( 2 _ r j J + (7r7 + 29)Li 2 (77) - -7t 2 (29t7 + 7) 



+ 18 log(2) log(2 - 77) + 2(47/ + 5) log(l - 77) log(77) - 18 log(2) logfa) 

1 /«n 2 ,0 _\ , / >. (2O77 3 + 207777 2 - 622177 + 4478) (1 - 77) 
+- (9O77 2 - 4877 + 17) log(77) + ^-L _ -L IS '1 

3 (33t7 3 - 11377 2 + IO677 + 4) (1 - r]f log(l - 77) 

(2 - 77) 2 77 ' ( ] 
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and for 1 P^\ 



g 1 [ 1 Pl 8] ] = 16(1 - V ) 2 ( 21og(l - V ) - log 



--(12^-6^ + 1) logfa), 



[i p [ 8 ] 
ri dM 



^pif 1 ] = 0, 



g 3 [ L PP] = 10(1 -vV 21og(l-7/)-lo, 



4m;? 



Ma 
4m 2 . 



- (8?7 2 - 85r/ + 119) (1 - 77) 



(29r? 2 - 244r/ + 347) (1 - 77) 



-- (30r/ 2 - 1577 + 7) logfa). 
o 



(62) 
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